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TR-86-6 

ERROR  PROBABILITY  FOR  DIRECT-SEQUENCE  SPREAD-SPECTRUM 
MULTIPLE-ACCESS  COMMUNICATIONS  OVER  NONSELECTIVE  AND 
FREQUENCY-SELECTIVE  RICIAN  FADING  CHANNELS 

Evaggelos  Geraniotis,  Member  IEEE 

Electrical  Engineering  Department 
and  Systems  Research  Center 
University  of  Maryland 
College  Park,  MD  20742 

ABSTRACT 

An  accurate  approximation  is  obtained  for  the  average  probability  of  error  in  an 
asynchronous  binary  direct-sequence  spread-spectrum  multiple-  access  communications 
system  operating  over  nonselective  and  frequency-  selective  Rician  fading  channels.  The 
approximation  is  based  on  the  integration  of  the  characteristic  function  of  the  multiple- 
access  interference  which  now  consists  of  specular  and  scatter  components.  For  non¬ 
selective  fading  the  amount  of  computation  required  to  evaluate  this  approximation 
grows  linearly  with  the  product  KN,  where  K  is  the  number  of  simultaneous 

transmitters  and  N  is  the  number  of  chips  per  bit.  For  frequency-selective  fading  the 

2 

computational  effort  grows  linearly  with  the  product  KN  .  The  resulting  probability  of 
error  is  also  compared  with  an  approximation  based  on  the  signal- to- noise  ratio.  Numer¬ 
ical  results  are  presented  for  specific  chip  waveforms  and  signature  sequences. 


This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under  contract  N00014-84-K-0023  and  in  part 
by  the  Systems  Research  Center  at  the  University  of  Maryland,  College  Park,  through  National  Science  Foundation 
CDR-85-00108. 


I.  INTRODUCTION 


This  paper  is  concerned  with  the  average  probability  of  error  for 
asynchronous  binary  direct-sequence  spread-spectrum  multiple-access 
(DS/SSMA)  communications  over  nonselective  and  frequency-selective  Rician 
fading  channels.  The  results  presented  represent  a  generalization  of  the 
analysis  of  [7]  which  considered  only  additive  white  Gaussian  noise  (AWGN) 
channels . 

The  system  model  is  described  in  detail  in  [6].  Throughout  this  paper 
we  restrict  attention  to  binary  DS/SSMA  systems  with  phase-shift-keyed  (PSK) 
modulation  and  chip  waveforms  which  are  pulses  of  arbitrary  shape  with 
duration  equal  to  the  inverse  of  the  chip  rate.  However  our  results  can  be 
extended  (using  the  results  of  [73)  to  include  other  forms  of  direct- 
sequence  modulation  like  quadriphase-shif t-keying  (QPSK),  offset  QPSK,  and 
minimum-shift-keying  (MSK). 

The  channel  model  considered  in  this  paper  is  the  Rician  or  specular- 
plus-Rayleigh  fading  channel  model  described  in  [1]  for  nonselective  fading, 
and  in  [5]  for  selective  fading.  In  the  second  case  the  fading  is  assumed 
to  be  frequency-selective  [33.  Although  our  results  can  be  extended  to 
time-selective,  as  well  as  doubly-selective  wide-sense-stationary 
uncorrelated-scattering  (WSSUS)  fading  channels  ([23,  [53),  we  restrict 
attention  to  frequency-selective  channels,  since  these  channels  often  arise 
in  practice  and,  as  shown  in  [53,  single  user  DS/spread-spectrum  systems 
outperform  biphase  PSK  systems  over  such  channels. 

In  [53  the  average  signal-to-noise  ratio  at  the  output  of  the  receiver 
of  a  binary  DS/SSMA  system  operating  over  time  and  frequency-selective 
Rician  fading  channels  was  evaluated.  In  [43  the  moment  space  bounding 


2 


technique  of  [12]  was  applied  to  bound  the  probability  of  error  of  such 
systems.  However  the  results  were  limited  because  of  unsurpassable 
difficulties  in  evaluating  high  order  moments  of  the  interference.  For  the 
same  reason,  the  approximation  based  on  the  Gauss  Quadrature  Rule,  which  was 
used  in  [10]  -  [11]  for  DS/SSMA  systems  over  AWGN  channels,  is  not 
computationally  attractive  in  this  case,  since  it  requires  calculation  of 
moments  which  involve  high  order  channel  autoeovariance  functions  that  may 
not  be  measureble  or  otherwise  available.  Similarly,  the  bounds  of  [133 
were  based  on  convexity  properties  of  the  complementary  error  function  which 
are  valid  only  in  the  AWGN  case.  Finally,  the  results  for  two  users  of  [14] 
for  nonselective  fading  channels  and  [15]  for  time-selective  fading  channels 
appear  difficult  to  extend  so  that  they  provide  an  efficient  and  accurate 
computational  method  for  the  multi-user  case. 

In  this  paper  we  obtain  an  approximation  with  a  computational 

requirement  which  for  nonselective  fading  is  linear  in  the  product  KN  and 

2 

for  frequency-selective  fading  is  linear  in  the  product  KN  ,  where  K  is  the 
number  of  simultaneous  transmitters  and  N  is  the  number  of  chips  per  bit. 
This  approximation  is  based  on  the  integration  of  the  characteristic 
function  of  the  multiple-access  interference  component  of  the  output  of  the 
correlation  receiver.  This  component  consists  now  of  a  fixed  (or  specular) 
part  and  a  random  (or  scatter)  part.  This  method,  which  we  refer  to  as  the 
characteristic-function  method,  was  applied  to  DS/SSMA  systems  in  [7]  and 
DS/SS  systems  with  specular  multipath  fading  in  [9]  where  it  was  shown  to 
provide  very  accurate  approximations  to  the  average  probability  of  error. 
Moreover,  as  shown  in  [8],  any  prespecified  degree  of  accuracy  can  be 


achieved  by  using  this  approximation  to  obtain  an  expansion  point  for  a 
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Taylor  series  representation  of  the  actual  probability  of  error.  In  the 
selective  fading  case,  however,  it  is  computationally  untractable  to 
evaluate  the  moments  of  the  interference  required  for  the  series  expansion 
method,  therefore,  in  the  sequel  we  consider  only  the  characteristic- 
function  method. 

We  should  also  point  out  at  this  point  that  the  DS/SSMA  system  under 
consideration  employs  conventional  matched  filter  receivers  combat  rather 
than  utilize  the  interference  (multiple-access  faded  and  nonf aded 
components)  by  discriminating  against  it.  This  is  in  contrast  to  the  work 
of  [17]  -  [18]  where  diversity  reception  is  used  to  combine  the 
contributions  of  the  faded  (reflected  in  those  cases) paths  and  use  them  in 
the  binary  decision. 

This  paper  is  organized  as  follows.  Nonselective  Rician  channels  are 
considered  in  Section  II,  and  frequency-selective  Rician  channels  are 
treated  in  Section  III.  In  each  of  these  sections  the  fading  channel  model 
is  first  presented.  Next  the  characteristic-function  method  is  applied  to 
each  particular  model  and  the  computational  requirements  are  considered. 
Finally  a  simpler  approximation  which  is  based  on  the  average  signal-to- 
noise  ratio  is  cited  for  the  sake  of  comparison.  In  Section  IV  numerical 


results  are  presented  for  specific  chip  waveforms  and  signature  sequences. 
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II.  PERFORMANCE  OF  DS/SSMA  COMMUNICATIONS  OVER  NONSELECTIVE  RICIAN  FADING 
CHANNELS 


A .  A  System  and  Channel  Model 

The  system  model  is  described  in  detail  in  [6]  and  [9].  Here  we  repeat  the 
basic  elements  of  this  model  so  that  the  concepts  and  notation  which  are 
necessary  for  the  rest  of  the  paper  are  introduced. 

The  k-th  transmitted  signal  for  a  binary  DS/SSMA  sy;stem  with  PSK 
modulation  and  arbitrary  chip  waveform  can  be  expressed  as 

s^Ct)  =  Re  (x^(t)  exp(  j  2-rrf^t)  }  (1) 

where  Re  stands  for  real  part,  j  =  /-l,  f  is  the  common  carrier  frequency, 
and  the  DS- spreaded  signal  x  (t)  is  given  by 

x^Ct)  =  /2P  bk(t)f(t)ak(t)  expCjO^).  (2) 


In  (2)  P  is  the  power  in  each  of  the  K  transmitted  signals;  the  equal  power 

assunption  is  made  for  convenience  in  presenting  numerical  results,  it  is  not 

necessary  for  the  methods  considered  in  this  paper.  The  phase  angle 

0^  introduced  by  the  modulator  is  uniformly  distributed  in  [0,2n  ].  The 

shaping  waveform  ^(t)  that  appears  in  (2)  is  periodic  with  period  T  (the 

c 

duration  of  a  chip),  it  is  defined  by  VCt)  =  ip(s)  for  s  =  t(mod  T  )  where  \p(s) 
is  a  chip  waveform  of  arbitrary  shape  (see  [6]  and  [9])  which  satisfies  a 
time-limiting  constraint  and  an  average  energy  constraint.  The  data  signal 
b  (t)  consists  of  rectangular  pulses  of  duration  T  which  take  values  +1  and  -1 

K. 


(k) 


The  information 


with  equal  probability.  The  £-th  pulse  has  amplitude  b 

A/ 

(k) 

sequence  (b  )  is  modeled  as  an  i.i.d.  sequence.  The  code  waveform  a  (t)  is 

/v  rC 

a  periodic  sequence  of  rectangular  pulses  of  amplitude  +1  or  -1  and  duration 


(k) 


We  assume  that  there  are  N  code 


T^.  The  j-th  code  pulse  has  amplitude  a 

pulses  in  each  data  pulse  (T  =  NT^)  and  that  the  period  of  the  signature 
(k) 

sequence  (a,  )  is  N. 

1 
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For  a  system  with  K  asynchronous  simultaneous  transmitted  signals^  the 

i/jjercsl’ed 

received  signal  at  a  receiver  which  ispa!  the  i-th  signal  is  given  by 
K 

r(t)  =  I  y,(t-T  )  +  n(^t^  O) 

k-1 

In  (3)  is  the  time  delay  for  the  communication  link  between  the  k-th 
transmitter  and  the  i-th  receiver.  The  process  n(t)  is  a  white  Gaussian  noise 
process  with  two-sided  spectral  density  N  /2.  Finally,  y  (t)  is  the  fading 

U  K. 

channel  output  to  input  s^(t). 

A  Rician  nonselective  fading  channel  is  described  by  the  following 
input-output  relationship 

yfc(t)  =  s^t)  +  Re  {\A£k^  exp(  j0*^  ^(t)  exp(  j2irf^t) }  (4) 

for  SLT  <  t  <  (X,  +  1)T.  In  (4)  y^  is  a  nonnegative  real  number  and  the 

(k) 

nonnegative  random  variable  A  satisfies  the  normalization  constraint 

X/ 

E{  [A^k)  ]2}  =  1.  (5) 

The  attenuation  of  the  signal  strength  due  to  the  fading  during  the  time 

(k) 

interval  [JIT,  (£  +1)T)  is  thus  represented  by  y.  Av  ,  and  the  phase  shift  due  to 

K.  X/ 

(k) 

the  fading  is  denoted  by  0  .  Therefore,  for  y  >  0  the  output  signal  is  the 

Xj  K. 

sum  of  a  nonfaded  version  of  the  input  signal  (specular  component)  and  a 

(k) 

nondelayed  faded  version  of  the  input  signal  (scatter  component).  For 

(k)  r  i 

Rayleigh  distributed  and  0^  uniformly  distributed  in  [0,2ir  j  this  model  is 
discussed  in  [l].  Notice  that  from  (4)  and  (3)  the  channel  model  reduces  to  the 
AWGN  model  if  y  =  0.  Similarly  for  y  ■*  “  (very  weak  specular  component  in 

K.  K. 

comparison  to  the  faded  component)  the  channel  reduces  to  the  Rayleigh 
nonselective  channel. 

Regarding  the  statistical  dependence  of  the  attenuation  in  different  data 

(k)  (k) 

bit  intervals,  we  will  consider  the  cases:  (i)  A„  =  A 
’  %  m 


for  several  values 
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of  adjacent  %  and  m  (the  length  of  the  burst  is  immaterial  for  our  analysis  as 

(k)  (k) 

long  as  it  is  longer  than  three  adjacent  bits)  and  (ii)  A^  and  are 

independent  if  H  *  m.  Case  (i)  corresponds  to  a  system  with  no  interleaving 

(or  partial  interleaving)  and  a  channel  with  fading  statistics  which  remain 

invariant  over  the  duration  of  several  data  bits.  An  example  of  case  (ii) 

arises  in  a  system  which  is  fully-interleaved  or  a  fast  fading  channel.  Similar 

(k) 

restrictions  are  imposed  on  the  phase  sequence  (0^  ). 

The  i-th  receiver  is  matched  to  the  i-th  transmitter  signal;  it  is  assumed 
capable  of  acquiring  time  and  phase  synchronization  (coherent  demodulation)  with 
the  nonfaded  component  of  the  i-th  signal.  Therefore,  we  may  measure  time 
delays  and  phase  angles  relative  to  and  6^,  respectively,  and  thus  set 

=  0^  =  0.  If  there  is  not  a  nonfaded  component  we  assume  that  the  phase  of 
the  faded  component  of  the  Rician  nonselective  fading  channel  [see  (3)  ]  changes 
very  slowly  so  that  the  receiver  can  again  acquire  time  and  phase  synchroniza¬ 
tion.  In  this  case  we  should  let  0.  +  gJ*^  =  0. 

l  0 

Once  time  and  phase  acquisition  is  completed  the  received  signal  is  passed 
through  a  correlation  receiver  which  outputs 


Z.  =  /  r(t) 'F(t)a  ,(t)  cos(2xf  t)dt. 
1  0  1 


(6) 


If  Z  is  positive  the  receiver  will  decide  that  a  1  was  sent,  otherwise  it  will 
output  a  -1.  In  practical  s pread-s pectrun  systems,  the  common  carrier  frequency 
f  and  bit  duration  T  are  such  that  f ^  >>  T  consequently,  the  double-frequency 
terms  in  the  integrand  of  (6)  may  be  ignored. 

Finally,  concerning  the  interaction  among  the  different  communication  links 
between  the  users,  we  may  measure  the  delays  T^f  espect  to  ,  set  =  0,  and 
assume  that  for  k  A  i  (mod  T)  is  uniformly  distributed  in  [0,T].  Then  the 
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(k) 

delays  t,  ,  the  phase  angles  0  ,  and  the  data  streams  (b^  )  are  assumed  to  be 

K  K 

mutually  independent  for  a  given  k ,  as  well  as  for  any  two  different  transmit¬ 
ters.  Finally,  we  assume  that  all  communication  links  fade  independently;  there¬ 
fore,  all  the  random  variables  and/or  random  processes  which  characterize  the 
fading  are  independent  for  any  two  different  links. 


B.  Evaluation  of  Error  Probability  via  the  Characteristic-Function 
Method 

By  ignoring  the  double-frequency  terms  in  (6)  we  may  write  as 
Z.  =  (1/2  EbT)1/2  (b^X)  +  +  /.),  (7) 

In  (7)  =  PT  is  the  received  energy  per  bit  for  the  nonfaded  component, 

the  first  term  inside  the  parenthesis  is  the  desired  signal  component  (i.e., 
the  i-th  information  signal  for  0  <_  t  <  T),  represents  the  normalized 
contribution  of  the  AWGN  and  of  the  faded  component  of  the  i-th  signal,  and 
ij,  stands  for  the  normalized  multiple-access  interference  (which  now  has 
nonfaded  and  faded  components)  due  to  all  the  other  users  besides  the  i-th 
user.  We  actually  have  that 


n.  =  n  +  Y .  F . 
i  ii 


(8) 


and 


h  ■  i  (Ik,i  *  Vm1’ 

kri 


(9) 


where  n  is  a  zero-mean  Gaussian  random  variable  with  variance  NQ/2Eb,  is 

the  normalized  interference  due  to  the  faded  component  of  the  i-th  signal, 

I,  .  and  F,  .  represent  normalized  interference  due  to  the  nonfaded  and 
k ,  l  k  ,  i 

faded  components  of  the  k-th  (k/i)  signal. 

Due  to  the  additive  form  of  the  multiple-access  interference  [see  (9)], 
the  c’naracteristic-f unction  method  is  recommended  for  the  evaluation  of  the 
probability  of  error  for  the  DS/SSMA  systems  under  consideration.  As  in 
[7]  we  can  write  P  .,  the  average  probability  of  error  for  the  i-th 

0  a  1 
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receiver  (the  averge  to  be  considered  with  respect  to  all  te  random 

variables  introduced  by  the  transmitters  and  the  channel)  as 

P  .  =1/2  P { Z .  <  0|b^l)  =  +1}  +  1/2P{Z.  >  0|b^l)  =  -1} 
e,i  l  -  1  0  l  1  0 

-  1/2  -  i/2P { — i  <  n.  +  /.  <  i }. 

where  p.  and  1 .  were  defined  above.  Let  d>  and  6.  denote  the 

l  i  rn-  i 


characteristic  functions  of  the  random  variables  of  and  i.  ,  respectively. 

Since  the  distributions  of  and  considered  in  this  paper  are 

symmetric,  the  characteristic  functions  <j>  and  <|> .  are  real-valued  even 

qi 

functions.  Therefore  as  in  [7]  -  [9]  we  can  write 
“  -0  "If00  “I,  .  .  .  ,  .  r.  .  /  nt, 


.  =  P  .  +  tt  u  (sinu)  <j>  (u)  (u)  ]du 

,i  e,i  }  n,  i 


where  P  .  is  the  average  probability  of  error  in  the  absence  of  multiple- 
e ,  l 

access  interference  (i.e.,  when  [.  =  0  or  K  =  1).  Because  of  (9)  and  the 

independence  assumptions  of  Section  II. A  above  we  can  write  for  each  real  u 

4,.(u)  =  n  E{exp[ju(Ik  .  +  YRF  .)]},  (12) 

k*i  ’  ’ 

which  facilitates  the  computation  of  the  characteristic  function  of  the 
multiple-access  interference. 

For  the  Rician  nonselective  fading  channel  of  Section  II. A  we  can 


easily  derive  that 

F.  -  b^A^coMO*1'), 

Vi  -  T"1[b(,k)Vi(v  *  t4k)Vi(v:|o°s8k' 


(13a) 

(13b) 


Vi  -  T'1[t’ik)#lk)  Vi(V003<8k 8  9,(k>) 

*  b2k)4k)  “k,i(\)ob3(9k  *  4k))]' 


(13c) 


where 


3.'  -  8  -  2irf  t,  . 

k  k  c  k 
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Recall  that  x .  =  6 .  =  0  so  that  0! 

li  l 


0.  The  arguments  of  the  partial 


continuous  crosscorrelation  functions  R  .  and  R  .  (defined  in  [6])  should 

K  |  1  K  f  1 

(k)  (k) 

be  considered  modulo  T.  Let  b  =  (b  ,  b  )  denote  the  pair  of  the 

K  I 

(X^-l)-th  and  Xk~th  consecutive  data  bits  (suppose  that  for 

(k) 

k  /  i  XT  <  x  <  (X  +  1 ) T) .  If,  for  simplicity,  we  assume  that  the  A  's 

K  K  K  X» 


are  Rayleigh  distributed  in  [0,  2ir],  then  since  6^  +  0^ 


(k) 


(mod  2tt)  is  also 


uniformly  distributed  in  [0,  2tt]  the  random  variables  F.  and  F  .  are 

1  K  f  1 

conditionally  zero-mean  Gaussian  when  conditioned  on  b^^  and  (bk,  x^), 


respectively. 

Next  we  proceed  to  the  evaluation  of  the  characteristic  functions  of 
and  which  are  necessary  for  the  computation  in  (11).  Let  a2  denote  the 


variance  of  defined  by  (8) ,  then 

4  ■  (2Eb/Nor'  * i/2y!’ 

pe.l  *  Q(h>. 

where  Q  is  the  complementary  error  function,  and 

2  2 

d>  (u)  =  exp  (-1  /2u  a . ) . 
rni  i 


(14) 

(15) 

(16) 


To  proceed  further  we  need  to  evaluate  the  variance  of  the  conditionally 

ssian  random  variable  F  .  given  by  (13c).  It  turns 

k ,  i 

out  that  for  (i)  constant  fading 


Gau 


',ar(Fk,i|W  ■  i/2  'I'“2[»Sk>Rk,i(V  * 


and  for  (ii)  independent  fading 
Var  {F 


,  .  lx,  }  =  1/2T  2{ Rf  .(t.  )  +  Rf  .(t,  )}- 

k.i1  k  k,i  k  k,i  k 


For  £T  <  x  <  (£+1  )T  and  any  arbitrary  function  h  define 
c  —  c 


a2(£,h,x)  =  1/2  T  2[h(£+1 )R  (t)  +  h(£)R^(x)]2, 

of  .  ( £ ;  x )  =  1/2  T-W  (x)  +  R2  (t)], 

K  ,  1  K  ,  1  K  ,  1 


(17) 

(18) 

(19) 

(20) 


and 


10 


f  (u;  ft ,h;  t )  =  |  |7r/2  cos{^[h(S,+1  )R^(x)  + 


h(J,)R^(T)]cos0}d0.  (21) 


In  (19)  and  (21)  R.  and  R,  are  the  continuous  partial  autocorrelation 
'Ij  ^ 

functions  of  the  chip  waveform.  The  functions  R  .  and  R  .  depend  on  them 

K  j  1  K  |  1 

and  on  the  discrete  crosscorrelation  function  of  the  k-th  and  i-th  signature 
sequences  C  . .  Then  after  some  straightforward  but  cumbersome 

K  ,  1 

manipulations,  (12)  becomes 

N-1  T 

Mu)  =  H  {(2N)"1  l  T"1  f  °  {f(u;l,0.  .  ;T)exp[-l/2u2Y2a2()l,e|  .;x)] 
k^i  £=0  c  J  0  k>1  k  k>1 

+  f(u;lt0k  .;T)exp[-1/2u2Y2a2U,ek  .;T)}di} 

(22a) 


for  (i)  constant  fading,  and 

N-1  T 

Mu)  =  n  {  ( 2N) _  1  I  t”1  f  °  [f(u;fc,0  ;t)  +  f(u;l,e.  .;t)] 

1  k^i  S>=0  J0  k’x 

•exp[-l/2u2Y2a2  ± ( S. ;  t ) ]dx }  (22b) 

for  (ii)  independent  fading. 

If  the  channel  encountered  by  the  i-th  transmitted  signal  is  a  Rayleigh 

nonselective  fading  channel  (i.e,  there  is  no  non-faded  component  of  the 

i-th  signal)  we  need  to  modify  the  above  results  in  the  following  way.  We 

start  with  finite  Y.  and  define  E,  =  (1  +  Y2)PT  to  be  the  total  received 

i  b  l 

energy  at  the  i-th  receiver  in  the  absence  of  AWGN  and  multiple-access 

-2  2  2 

interference  (i.e.,  K  =  1),  and  Y  =  (1  +  Y  )/(1  +  Y.)  to  be  the  ratio  of 

K  K  1 

the  total  received  energy  from  the  k-th  signal  to  the  total  received  energy 
from  the  i-th  signal.  For  large  Y.  ,  we  assume  (see  Section  II. A)  that  the 
receiver  acquires  time  and  phase  synchronization  with  the  faded  component  of 
the  i-th  signal,  so  that  we  may  set  9^  +  0 ^  =  0  and  together  with  the 
definitions  above  write  (7)  and  (8)  -  (9)  as 


11 


1  /? 

Z.  =  (1/2E.T)  [n+ 

i  b 


/l 


+Y . 


(b<1)=ose1+'r1b'i)A<i,)+I  - 
2  W  / 


k 


+Y. 


(Ik,l*YkFk,i)]' 


(23) 


If  we  fix  E,  and  let  Y.  ■*  00  (23)  reduces  to 
b  l 


Z.  =  (l/2EbT) 
(i) 


(24) 


where  is  Rayleigh  distributed  (with  second  order  moment  1),  =  n  is 

zero-mean  Gaussian  with  variance  N./2E  and  I .  is  given  by 

Ob  i 


£  -  l 


Y, 


k/i 


/l  +Y2 


(I,  .  +  Y  F  .) 
k ,  x  k  k ,  x 


(25) 


where  I,  .  and  F.  .  are  as  defined  in  (13b)  -  (13c). 

K  ,  1  rC  f  X 

A  discussion  of  the  computational  requirements  for  formulas  of  the  type 
of  (11)  can  be  found  in  [7]  -  [8]  and  since  the  situation  here  is  very 
similar  with  that  treated  there,  it  will  not  be  repeated.  However,  it  is 
worth  observing  that  since  there  is  a  summation  of  N  terms  and  a  product  of 
K  terms  in  (22a)  -  (22b)  the  computational  effort  for  this  method  grows 
linearly  with  the  product  KN. 

Starting  from  (24)  and  proceeding  as  for  (10),  we  find  that  P  (A  J"  ) 
the  conditional  probability  of  error  given  A^^  can  be  expressed  initially 
as 

Pg  .(A^)  =  1/2  -  1/2  P{-Aq1}  <  n.  +  -T.  <  Aq1)}’  (2&) 

and  finally  [cf.  (11)]  as 

P  .U{ni])  =  P°  .(A^)  +  it"1  f  u_1sin(A^l)u)(!)  (u)[1  -  <b.(u)]du. 

e,i  0  e ,  i  0  J0  0  ni  i 

(27) 


Since  n.  is  now  zero-mean  Gaussian  with  variance 
i 

a2  =  (2E./Nn)_1 
l  b  0 


(28) 
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2 

4>  (u)  is  given  by  (16)  with  o.  defined  by  (28)  and 

^i 

pe,i(A01>)  ‘  Q<A01)/Oi>-  <29) 

To  obtain  the  average  probability  of  error  P  ,  we  should  average  the  second 

e ,  l 

member  of  (27)  with  respect  to  the  Rayleigh  distributed  random  variable 

The  result  (see  Appendix  A  of  [9]  for  the  proof)  is 

P  .  =  P°  .  +  1/2tt_1/2  f“  exp  (-1/Mu2  )4»  (u)[1  -  4> .  ( u )  ]du  (30) 

e,i  e ,  i  JQ  ^  i 

where 

P°  .  =  1 / 2 { 1  —  [ 1  +  (Eb/NQ)_1 ]_1/2}  (31) 

and  (u)  is  given  by  (22a)  and  (23b)  for  case  (i)  constant  fading  and  (ii) 

2 

independent  fading,  respectively,  evaluated  [because  of  (25)]  at  Yku//l+Yk 
instead  of  u. 
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In  (34)  the  function  p  . (n)  (see  [6])  depends  on  the  discrete 

k ,  l 

crosscorrelation  function  C,  .  while  m,  and  m'  (again  see  [6])  depend  on  R 

k,i  ip  ip  >J 

and  R  .  For  Rayleigh  nonselective  fading  (33)  becomes 

SNR.  -  «VV-’  *  j.  ^  T-3[„kji(0)V„kjl(,)m;]r1/2  (36) 

and  the  resulting  approximation  is 

P°  .  =  1/2{1-C1+(SNR.)_2]_1/2}.  (37) 

Notice  that,  in  contrast  to  the  approximation  based  on  the  characterisitic- 
f unction  method,  the  approximations  obtained  via  the  signal-to-noise  ratio 


do  not  ditinguish  between  the  constant  and  independent  fading  cases 


considered  in  Section  II. A. 
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III.  PERFORMANCE  OF  DS/SS  COMMUNICATIONS  OVER  FREQUENCY-SELECTIVE  RICIAN 
FADING  CHANNELS 

A .  Channel  Model 

The  input-output  relationship  for  a  wide-sense-stationary  uncorrelated- 
scattering  (WSSUS)  frequency-selective  fading  channel  can  be  expressed  as 
yk(t)  =  Sk(t)  +  Re^uk(t-^kT)exPU2Trfc(t-vkT)]}  (38) 

where 

Uk*'t)  =  \  f  hk^t>Oxk(t-C)dC-  (39) 

j  —CO 

The  fading  process  h  (t,c)  (which  can  be  thought  of  as  the  time-varying 
impulse  response  of  a  low  pass  filter)  is  a  zero-mean  complex  Gaussian 
random  process  wich  has  autocovariance 

E{hk(t,s)h*(s,£)}  =  gk(0<SU-£),  (40a) 

and  satisfies  the  normalization  constraint 

f  gR(c)dC  =  1,  (40b) 

'  -CO 

and  the  necessary  condition  for  stationary  bandpass  processes  (see  [2],  and 
the  discussion  in  Section  2.2  of  [4]) 

E{hk(t,c)hk(s,5)}  =  0.  (40c) 

The  first  term  of  (38)  is  again  termed  the  specular  component  while  the 
second  term  in  (38)  is  now  called  the  diffuse  faded  component.  If  Yk  is  a 

positive  real  number,  the  channel  is  a  Rician  fading  channel.  If  Y  ■+  <*> 

K 

the  faded  component  becomes  dominant  and  the  channel  is  a  Rayleigh  fading 
channel.  This  model  is  basically  described  in  detail  in  [2]  -  [3].  in 
particular  the  Rician  fading  channel  is  described  in  [4]  -  [5].  Here  we 
consider  a  slightly  more  general  model  since  for  the  Rician  case  we  allow 
[see  (38)]  an  average  delay  vkT  (vk  is  a  nonnegative  integer)  between  the 
nonfaded  specular  component  and  the  diffused  faded  component. 
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The  WSSUS  fading  model  can  be  thought  of  arising  when  the  transmitted 
signal  encounters  a  slowly-moving  scattering  random  medium  which  can  be 
modeled  as  a  layered  scatterer  consisting  of  a  large  number  of  layers  of 
infinitesimal  thickness.  It,  is  therefore,  a  frequency-selective  channel, 
i.e.,  it  is  dispersive  only  in  frequency.  The  high-frequency  (HF) 
ionospheric  and  the  microwave  tropospheric  scatter  channels  fit  the  above 
description;  Section  2.1  of  [4]  provides  more  examples  of  such  channels. 

We  impose  a  limitation  on  the  selectivity  of  the  channel  in  order  to 
facilitate  subsequent  analysis.  In  particular,  as  in  [5]  and  [3]  it  is 
assumed  that 

gRU)  =  0  for  |c|  >  T  (41) 

which  is  a  constraint  on  the  frequency  selectivity  of  the  channel  that 
allows  us  to  restrict  attention  to  the  intersymbol  interference  from  the  two 
adjacent  data  bits. 

With  regard  to  the  average  delay  v^T  between  the  nonfaded  specular 
component  and  the  faded  scatter  component  of  yR(t)  in  (38),  we  will  consider 
the  two  extreme  cases  (i)  vR  =  0  (no  delay)  and  (ii)  vR  ^  3  (large  delay). 
The  choice  of  the  integer  3  as  the  lower  bound  for  v.  in  (ii)  is  justified 

K 

by  the  fact  that,  in  view  of  ( 4 1 )  ,  it  uncorrelates  the  contributions  of  the 
nonfaded  and  faded  components  at  the  output  of  the  receiver. 

B .  Evaluation  of  Error  Probability  via  the  Characteristic-Function 
Method 

The  output  of  the  matched  filter  Z.  is  again  given  by  (7)  where  and 

are  given  by  (8)  and  (9),  respectively.  The  random  variable  q  in  (8)  is 

again  zero-mean  Gaussian  with  variance  NQ/2Eb;  the  nonfaded  component  of  the 

interference  (due  to  the  k-th  signal)  I,  .  is  given  by  (13b),  but  the  faded 

k ,  1 
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components  F.  and  F  .  of  (8)  -  (9)  are  now  found  [combine  (38)  -  (39)  with, 

1  K  j  1 

(3)  and  (6)]  to  be 


Re{F.} , 

T_1  f°°  h. (t-v.T1?)ri(t)?+v.T)exp(jOi)dcdt, 

*  — 03 


(t"VvkT^)rk,i 


(t ,c+Tk+vkT)exp(jOk)dcdt , 


(42a) 

(42b) 

(42c) 

(42d) 


where 


r,  .(t,c)  -  a,  (t-c)a.(t)b.  (t-c)b.(t)¥(t-c)'P(t) 
k ,  i  k  i  k  x 

and 


$k  =  9k  ’  2-fc(W)* 

The  function  r  .  is  denoted  by  r. .  Since  h  is  a  complex  Gaussian  random 

1,1  i  x 

field,  the  random  variables  F.  and  F.  .  are  conditionally  Gaussian  random 

1  K  f  1 

variables  (conditioned  on  time  delays,  phase  angles,  and  data  streams). 

Recall  that  v  T  (1  <  k  <  K)  accounts  for  the  average  delay  between  the 
k 

nonfaded  and  the  diffused  components  of  the  channel  output  (see  Section 


III. A) 

Next,  we  proceed  to  the  evaluation  of  the  average  error  probability 
p  via  the  characteristic-function  method.  As  in  the  nonselective  fading 

e ,  i 

case,  we  need  to  evaluate  Pp  ,  (j)  ,  and  <j>.  for  use  in  (11).  We  follow  the 

same  steps  as  before.  The  random  variable  F ^  is  zero-mean  conditionally 
Gaussian  with  conditional  variance 
Var{F.|b!}  =  1/2  E(F.F*} 

=  T-2  [T  g.(e)[  R2U)  +  Rj(c)  +  (b1(J)+b2t))RjL(?)R1U)]dC.  (43) 

Jo  1 

where  b!  =  (b^k  b^k  denotes  the  pair  of  the  (v^-1  )-th  and  v^-th 
consecutive  data  bits  of  the  i-th  data  stream.  To  derive  (43)  we  used 
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(42a),  (42b),  (40a)  -  (40c),  and  (41).  The  Gaussian  random  variable  n  has 
then  conditional  variance 

V ar  { n  .  I  b }  =  (2EK/N  )-1  +  Y2Var  {F  .  I  b } .  (44) 

i 1  — i  b  0  l  i'-i 


Then,  if  we  define 

[a^]2  =  (2E./N)"1  +  Y2  fTg.  (?)[R2(c)+R2(^)]d?  (45a) 

l  b  0  1  J0  1  1 

[o[2)]2=  (2Eb/NQ)_1  +  T2T_2|T  gi(?)[Ri(c)+Ri(c)l2dc  (45b) 

[af3)]2  =  (2Eb/NQ)_1  +  Y2T“2JT  gi(c)[Ri(4)-Ri(4)]2dC,  (45c) 

we  can  write 


P 


n 

e 


=  I  cQCCc^]"1) 

j  =  1  J 


(46) 


and 

3 

<f>  =  l  c  exp{-1/2[ua[J)]2}  (47) 

ai  j=1  J 

where  c^  =1/2  and  c2  =  c^  =  1/4.  Note  that  these  two  results  are 

independent  of  the  events  =  0  (no  delay:  case  (i)  of  Section  III. A  or 

v.  >  3  (large  delay:  case  (ii)  of  Section  III. A. 

The  next  step  is  to  evaluate  <j> .  through  (12)  with  I.  .  and  F.  .  given 

1  K  j  1  K  j  1 

by  (13b)  and  (42c)  -  (42d),  respectively.  We  first  compute  the  conditional 

variance  of  F,  . ,  this  turns  out  to  be 

k,i  ’ 

Var{Fk,il-k’T}  =  1/2T"2{|X  gk(4-T+T)[b^)Rkj.(0+b^)Rk>.(i;)]2dc 

iTg  8k(t”T)[b2')RK,iU)  *  b3'>RK,iU)]2dC 
+  jT  gk(^-x-T)[b^)Rk  .(c)  +  b^}Rk  i(c)]2d?}. 


(48) 

In  (48)  =  (b^,^,  b^^,  b  b^,1',<^)  (k  *  i)  is  the  quadruple  which 

consists  of  the  (p  -2)-th,  (p.-l)-th,  p  -th,  and  (p  +1 )-th  consecutive  bits 

K  K  K  K 

of  the  k-th  data  stream,  pk  =  Xk  +  vk,  and  x  =  xk(mod  T)  where 
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£T  <  x  <  (£+1 )T  .  The  conditional  variance  about  takes  on  8  different 
c  —  c 

values  which  we  denote  by  [a.^?  (£;x)]2  and  [af^  (£;x)]2  for  j  =  1,2, 3, 4. 

K  f  1  K  f  1 

Then  the  characteristic  function  <j> .  is  derived  as 

vi 


i.  (u)  =  n  { (8N) 
k*i 


-1 


N-1 


l  T1  C[f(u;£,e.  .  ;x)  l  expf-1/2[uTVj;U;i 

K,I  KK,1 


(j), 


)rj 


£=0 


j  =  1 


+ 


4 

f (u; £,  8  .  ;  x )  l 

>1 


exp{-1/2[uYko^(£;T)]2}]dT} 

(49a) 


for  (i)  no  delay  (v 

<t>.(u)  =  n  { (8N)_1 
1  k*i 


=  0 ,  all  k) ,  and 

N-1  _  T 

l  T  f  °  [f (u; £ ,  0 
£=0  J0 


+  f (u;£,0  . ; t )  ] 

K  ,  1 


•  I  [exp {-1  /2[uY  of.^!*(£;T)]2} 
j-1  k  k’ 

+  exp{-1/2[uYRak^  (£;t)]2}]dx}  (49b) 
for  (ii)  large  delay  (v,  >  3,  all  k).  In  case  (i)  I.  .  [given  by  (13b)] 

K  —  K  j  1 

and  F  .  above  are  correlated  [  t> '  =  (bv  ,  b  ,  b,, ,  )].  In  case  (ii)  I  . 

Kjl  KIKH  K , 1 

and  F  .  are  uncorrelated  (b  and  b'  have  no  components  in  common).  These 

K  f  1  K  K 

facts  are  reflected  in  the  form  of  <{k  in  (49a)  and  (49b)  (the  function  f  and 
the  sum  with  the  four  terms  are  the  characteristic  functions  of  I  .  and 

K  ,  1 

F  .  given  x  respectively).  The  function  f  is  as  defined  in  (21). 
k ,  i 

2 

To  see  on  the  computational  effort  grows  linearly  with  KN  ,  notice  that 

in  (49a)  -  (49b),  besides  the  product  of  K-1  terms  and  the  summation  of  N 

terms,  the  integrals  split  into  summations  with  a  total  of  N  terms  as  well 
£=1  (n+1 )T 

[e.g.,  (T  =  I  f  C  +  [T  ,  for  £T  <  x  <  U+1  )T  ]. 

J0  n=0  JnT  J  £T 

c  c 

C.  Approximation  Based  on  the  Average  Signal-to-Noise  Ratio 


The  average  signal-to-noise  ratio  in  this  case  is  given  by  [cf.  (34)3 
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SNR.  =  {  [  a  f 1  }  ] 2  +  I  (m2)[y.  •  ( 0 )  m  +  y  .(1)m']}  1/2  (50) 

1  1  ^  K  K,1  K ,  1  IjJ 

( i ) 

where  ck  is  given  by  (45a)  above.  This  is  a  slight  generalization  of  the 


results  of  [5]  for  arbitrary  chip  waveforms.  Notice  that  SNR^  and  the 
resulting  approximation  (35)  cannot  distinguish  between  cases  (i)  and  (ii) 


(no  delay  versus  large  delay)  of  the  fading  model  of  Section  III. A. 
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IV.  NUMERICAL  RESULTS  AND  CONCLUSIONS 

All  the  signatures  sequences  employed  in  this  section  are  AO/LSE  m- 
sequences  of  lengths  N=31  and  127.  In  particular,  for  a  DS/SSMA  system  with 
K  users  the  first  K  sequences  of  Figures  A.1.(a)  ( N=  3 1 )  and  A.1.(b)  (N=127) 
of  [15]  are  used. 


For  P  .  ,the  approximation  to  the  error  probability  which  is  based  on 
e ,  x 

the  characteristic-function  method, we  use  Simson's  integration  technique 

-1  4 

with  the  same  parameters  as  in  Section  III  of  [7]  (L=20,  e<10  ,  n  =10, 

n  =10,  and  n=20) .  We  set  i=1  and  consider  P  for  the  error  probability  of 

T 

the  receiver  matched  to  the  first  signal. 

First, we  present  numerical  results  for  DS/SSMA  systems  with 
nonselective  Rician  and  Rayleigh  fading  channels.  In  Table  1  the  Gaussian 

_G  **”  ( i  ) 

approximation  to  the  error  probability  P  and  the  approximations  P  and 

e ,  i  q  ,  i 

obtained  by  the  characteristic-function  method  are  compared  for  a 

e ,  i 

DS/SSMA  system  with  sequences  of  length  31  and  K  =  2,  3,  and  4  simultaneous 
users.  The  communications  link  between  the  k-th  transmitter  and  the 


receiver  matched  to  the  first  signal  is  a  Rician  nonselective  fading  channel 

2  2 

with  relative  power  of  the  faded  component  7^=7  =.1,  1<k<K.  The  error 
probabilities  P^j  and  correspond  to  the  cases  of  (i)  slow  fading  and 

(ii)  fast  fading  discussed  in  Section  II. A.  We  observe  that  the  Gaussian 


approximation  is  somewhat  optimistic  when  compared  with  the  more  accurate 
approximations  based  on  the  characteristic  function  but  still  in  qood 
agreement  with  them.  Notice  that  slow  fading  causes  a  slightly  worse 
performance  of  the  DS/SSMA  system  than  fast  fading.  Recall  that  the 
Gaussian  approximation  results  in  the  same  value  of  the  error  probability 
for  both  slow  and  fast  fading. 
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/  «  \ 

In  Figure  1  we  plotted  both  P  1  (for  slow  fading)  and  P  versus  E  /N 

G  j  I  G  ,  !  DU 

for  a  DS/SSMA  system  with  K=3,  N=31  ,  and  Rician  nonselective  fading.  The 

2 

relative  power  of  all  faded  components  are  Y  =.01,  .05,  .1,  .2,  and  .4. 

2 

Notice  that  as  Y  increases  the  system  degrades  gracefully.  For  a 
nonselective  fading  channels  DS/SS  modulation  enables  the  matched  filter 
receiver  to  discriminate  effectively  against  the  nonf aded  and  faded 
components  of  the  interfering  signals  but  not  against  the  faded  component  of 

the  desired  signal.  This  explains  why  the  error  probability  is  already 
-3  2 

larger  than  10  for  Y  =  .1,  .2,  and  .4.  Finally,  notice  that  the  accuracy 

2 

of  the  Gaussian  approximation  improves  as  Y  increases. 

In  Table  2  we  compare  the  performance  of  DS/SSMA  systems  with 
rectangular  and  sine  chip  waveforms.  The  sine  chip  waveform  used  is 

ip ( t )  =  /2  sin  (irt/T  )  ,  0  S  t  iS  T  .  The  system  parameters  are  K=3,  and 

G  C  _ 

N=31 .  The  Rician  nonselective  slow  fading  channel  parameter  takes  on  the 

2  2 
values  Y  =  .01  and  .1.  For  moderately  heavy  fading  (Y  =.1)  the  rectangular 

chip  waveform  results  in  a  slightly  better  performance  than  the  sine  chip 

_  Q 

waveform  as  both  the  Gaussian  approximation  P  and  the  more  accurate 

e ,  i 

approximation  indicate.  For  light  fading  (Y  =.01)  the  Gaussian 

e ,  i 

approximation  still  shows  better  performance  for  the  sine  chip  waveform  for 
all  values  of  Eb/NQ;  however,  the  more  accurate  approximation  indicates  that 
the  sine  chip  waveform  outperforms  the  rectangular  chip  waveform  only  for 
low  signal-to-noise  ratios;  the  situation  is  reversed  for  higher  signal-to- 
noise  ratios. 

In  Table  3  we  present  results  for  a  DS/SSMA  system  with  K=3,  N=3 1  and 

Rayleigh  nonselective  fading  channels  with  relative  power  of  faded 

-2 

components  taking  values  Y  =.01,  .1,  and  .2.  The  system  performance  has  now 
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degraded  considerably  and  higher  signal-to-noise  ratios  were  considered. 

Notice  that  the  Gaussian  approximation  is  conservative  and  for  large  values 
-2 

of  Y  is  off  the  actual  value  by  one  order  of  magnitude.  On  the  other  hand 
the  values  of  the  accurate  approximation  for  (i)  slow  fading  and  (ii)  fast 
fading  do  not  differ  in  the  first  four  significant  digits. 

Next,  we  consider  DS/SSMA  systems  operating  through  Rician  frequency- 
selective  channels.  In  Table  4  the  Gaussian  approximation  is  compared  to 
the  approximation  obtained  via  the  characteristic-function  method  for  the 
cases  (i)  zero  average  delay  between  nonfaded  and  faded  components  and  (ii) 
large  delay  (see  Section  III. A)  and  DS/SSMA  systems  with  K=3  and  N=31  and 

p 

127.  The  relative  power  of  the  diffuse  component  is  Y  =.4  for  all  signals. 
Notice  that  the  system  performance  for  case  (ii)gives  slightly  lower  error 
probability  than  case  (i)  due  to  the  more  extensive  randomization  involved. 
The  Gaussian  approximation  results  in  an  optimistic  estimate  of  the  error 
probability  for  all  cases  but  it  is  satisfactorily  close  to  the  more 
accurate  results.  Recall  that  the  Gaussian  approximation  takes  the  same 
value  for  the  af ormentioned  cases  (i)  and  (ii). 

In  Figure  2  and  are  plotted  versus  E,/Nn  for  DS/SSMA  systems 

G  ,  1  6, I  DU 

with  N= 3 1  and  K=2,  3.  and  4.  The  autocovariance  function  of  the  frequency- 

1  I  ^  I 

selective  WSSUS  channel  is  triangular:  g.  (t)  =  — (1  -  -L~)  for  |t|  £  T  and  0 

K  1  1 

2 

otherwise.  The  relative  power  of  all  the  faded  components  is  Y  =.4.  Notice 

-5 

that  for  an  error  probability  of  10  the  Gaussian  approximation  is  off  by 
almost  1  dB  for  K=3.  For  fixed  E^/N^  the  Gaussian  approximation  is 
satisfactory  for  all  cases  as  it  is  also  indicated  by  the  results  of  Table 


4. 


23 


Finally,  in  Figure  3  we  have  plotted  versus  E^/Nq  for  a  DS/SSMA 

system  with  K=3  and  N=31  and  relative  power  in  the  diffuse  component  taking 
2 

values  Y  =0.,  .1,  .2,  .3,  and  .4.  Notice  the  graceful  degradation  of  the 

2  2 

system  performance  as  Y  increases.  The  value  Y  =0.  corresponds  to  DS/SSMA 
with  an  AWGN  channel.  A  comparison  of  Figures  1  and  3  verifies  that  DS/SS 
modulation  is  much  more  efficient  against  frequency-selective  fading  than 


against  nonselective  fading. 
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Table  2.  Probability  of  error  for  DS/SSMA  systems  with  different  chip 
waveforms  and  Rician  nonselective  fading  (K=3,  N=  3 1 ) . 
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Table  3.  Probability  of  error  for  a  DS/SSMA  system  with  Rayleigh 
nonselective  fading  (K=3,  N=  3 1 ) . 


y2=. 

.01 

y2=. 

.1 

y2=. 

.2 

/No 

P(i> 

e,1 

p(ii ) 
e ,  1 

P(1> 

e ,  1 

p(ii) 

e,1 

?e,, 

-p'i) 
e ,  1 

p(ii) 
e ,  1 

10 

6.42 

6.43 

6.43 

(xIO-2) 

6.47 

6.51 

6.51 

(xl 0~2 ) 

6.53 

6 . 60 

6.60 

(x10~2) 

20 

2.35 

2.34 

2.34 

(xl 0~2 ) 

2.55 

2.43 

2.43 

(xIO-2) 

2.74 

2.53 

2.53 

(xIO"2) 

30 

0.85 

0.78 

0.78 

(x10~2) 

1.34 

0.87 

0.87 

(x10~2) 

1  .72 

0.98 

0.98 

(xIO-2) 

40 

4.37 

2.59 

2.59 

(xio"3) 

1.14 

0.35 

0.35 

(x10_2) 

1  .56 

0.46 

0.46 

(xIO-2) 

50 

3.70 

0.89 

0.89 

(x10~3) 

11.2 

1.83 

1.83 

(x10~3) 

15.6 

2.87 

2.87 

(xl 0~3 ) 

00 

36.2 

1  .07 

1  .07 

-4 

(xl  0  ) 

11.2 

1  .05 

1  .05 

(x10_3) 

15.6 

2.08 

2.09 

(xl  o'"3) 

Table  4.  Probability  of  error  for  a  DS/SSMA  system  with  frequency- 


selective 

f  ading 

(K=3,  Y2= 

.4). 

N= 

31 

N=1  27 

E  /N  ( dB) 
b  0 

—  r* 
pU  1 

e,  1 

e,  1 

p(i;> 

e ,  1 

P(1) 

e,1 

-p(n> 

e ,  1 

6 

6.34 

6.40 

6.40 

(xl 0_3 ) 

3-20 

3.20 

3.20 

(xIO  3) 

8 

1  .57 

1 .62 

1  .62 

(xIO-3) 

3.78 

3.82 

3.82 

-4 

(xIO 

10 

3.1 0 

3.42 

3.41 

-4 

(xIO  ) 

1  .83 

1  .97 

1.97 

(xIO”5) 

12 

5.54 

6.97 

6.95 

(x10_5) 

3-31 

3.89 

3.89 

(x10“7) 

14 

1  .05 

1  .61 

1  .59 

(xIO-5) 

1  .91 

3.22 

3.20 

(xl 0~9 ) 

1  6 

2.47 

4.75 

4.73 

(xIO-6) 

0.44 

1  .02 

1 .01 

(xl 0_1 ° ) 

Figure  Captions 


Figure  1.  Probability  of  error  for  a  DS/SSMA  system  with  Rician 

2 

nonselective  fading  (K=3,  N=31,  and  y  =.01,  .05,  .1,  .2,  and  .4) 

Figure  2.  Probability  of  error  for  a  DS/SSMA  system  with  frequency-selective 
fading  (N=31,  y^=.4,  and  K=2,3,  and  4). 

Figure  3.  Probability  of  error  for  a  DS/SSMA  system  with  frequency-selective 

2 

fading  (K=3 ,  N=31,  and  y  =0.,  .1,  .2,  .4,  and  .6) 
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Figure  1.  Probability  of  error  for  a  DS/SSMA  system  with  Rician 

nonselective  fading  (K  =  3,  N  =  31,  and  y2  =  .01,  .05,  .1,  .2,  and  .4) 


Figure  2.  Probability  of  error  for  a  DS/SSMA  system  with  frequency-selective  fading 
(N  =  31,  ®  .4, and  K  *»  2,3,  and  4). 
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Figure  3.  Probability  of  error  for  a  DS/SSMA  system  with  frequency-selective  fading 
(K  **  3,  N  «  31,  and  y^  “  0.,  .1,  .2,  .4,  and  .6) 


